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Abstract 

The coordinate of a harmonic oscihator is measured at a time 
chosen at random among three equiprobable instants: now, after one 
third of the period, or after two thirds. The (total) probabihty that 
the outcome is positive depends on the state of the oscihator. In the 
classical case the probability is ^ ± |, but in the quantum case it is 
0.50 ±0.21. 



Introduction 

The coordinate q{t) of a harmonic oscillator depends on the time t, 

lif) = % cos t + po sin t , 

where go^Po are the initial coordinate and momentum; the mass is assumed 
to equal 1, and the period — to equal 27r. In the classical setup go and po 
are numbers, while in the quantum setup they are operators, 

Q(t) = Qcost + Psint, [Q,P] = i 

(assuming also h = 1). Here is a question trivial in the classical setup but 
nontrivial in the quantum setup. We choose r at random from the three- 
element set {0, 27r/3, 47r/3} and check, whether g(r) > or not. The (total, 
unconditional) probability of the event g(r) > depends on the initial state 
of the oscillator. The question: what is the maximum of the probability over 
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all states? In the classical setup the probability is either 1/3 or 2/3. 
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In the quantum setup the probability is 

prob(^) = 1(^1^(0) + E(27r/3) + E{47i/3)\i:) ; 

here ijj is the state vector, E(t) = 9{Q(t)), and 9{q) = 1 if g > 0, otherwise 
0. Thus, we want to find the upper bound of the spectrum of the operator 
£^(0) + E{27t/3) + E{Att/3). The question is nontrivial, because the three 
terms do not commute. Note that we perform a quantum measurement only 
once (at a time chosen beforehand), therefore the impact of the measurement 
on the state is irrelevant. 

A numeric computation reported in Sect. [1] shows that the upper bound, 
sup^ prob ('?/'), of the spectrum is close to 0.71. A rigorous result of Sect. [2] 
states that sup^ prob (■?/') < 1. Sect. O generalizes this result to E{0) +E{s) + 



1 Using the Wigner quasi-distribution 

The Wigner function (or quasi-distribution density) : ^ M corre- 
sponding to a state vector ip has several equivalent (sometimes, up to a 
coefficient) definitions. The 'tomographic' definition (see also |i2j. Sect. 6.2, 
Th. 6.1), stipulating that returns correct one- dimensional distributions, 
is based on the equality 

(1.1) jj f{aq + bp)W4q,p)dqdp= {ij\f{aQ + bPm 

for all a, 6 G M and all bounded measurable functions / : M ^ M. However, 
there is a catch: need not be integrable, that is, // dq'dp = oo 

for some ip (fo^' example, ipi^) = {q"^ + l)~^^^e"'^ in the Schrodinger repre- 
sentation). Thus, even the most well-known relation JJ Wjij{q,p) dqdp = 1 
needs a careful interpretation! 
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One may interprete the left-hand side of fll.ip as 

hm J J f{aq + bp)W^{q,p)h{{bq — ap)6)dqdp 
for an appropriate weight function h such as 



h(x) = e ^ or h(x) 



1 — |x| for |x| < 1, 
otherwise. 



(It would be natural to put h{x) = 1 for |a;| < 1, otherwise 0. However, I do 
not know, whether this function fits or not.) Treated this way, f 1 1.11) holds 
for all ip. 

Alternatively, one may characterize the whole map ip ^ by two 
conditions: 

(a) there exists a dense set of state vectors ip such that is integrable 
and (11.11) is satisfied; 

(b) W^„{q,p) — ^ W^{q,p) (as n — >• oo) uniformly in q,p whenever ipn — i^- 
In the Schrodinger representation, 

W4q,p) = ^ y"^(g+|)^(g-|)e^^^da: 

for ip G -L2(M), and of course, 

iQi;){q) = qi;{q) , (PV^)(g) = -i^V^(g) 

for ip good enough. 
We have 



{ilj\E{t)\ij) = [[ e{qcost + psmt)W4q,p)dqdp= f^' W;''\^) d^ , 

J J Jt--K/2 



where 



oo 



W^^{(p) = / iy^(rcos(/9,r sinyj) rdr 







may be called the angular Wigner function. Thus, 
(1.2) 

prob(^) = ^(^|E(0) + E(27r/3) + ^(477/3)1^) = ^ f f {v)W;^\^) d^ , 



— TT 
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where /(■) takes on two values, 1 and 2, as follows: 

7r/2 

(1.3) 




The inequality | < prob(^/') < | holds if W^^{}^^ > for all (p, but fails in 
general. 

Eigenvectors of the Hamiltonian are the number states |0), |1), |2), . . . ; 
(Q2 + P2)|^) = (2n + l)\n). 

Generally, 



i> = '^Cn\n), ^|c„|^ = l, 

n=0 n=0 

{(1,P) , 



where functions Wm,n are defined by 

f{aq + bp)wm,n{(l^P) dgdp = {m\f{aQ + bP)\n) 

J J 

In fact, 

(1.4) Wrr.,n{q,p) = izll!!. H(2^)— e-2N'L(^— )(4|;^n 

TT V ^! 

for m < n (and Wn,m = Wm,n)', here z = (q+ip) / V^, and Lm^ is the associated 
Laguerre polynomial [H Sect. 10.12], 

iS'W = (-l)'5;;iWW, 



Equality (O) is a combination of P (3.30)] and [i (4.36)]. 



For the angular Wigner function, 



Wm,n{r COS (f, T siu if) r&T . 



Using (11 .4^ we get, first. 



and second, 
(1-5) <^.(0) 



Wm,n{r, 0) rdr 







^ / (y2r)"-"^e-"'L(;^-'")(2r')rdr 



-1 



vm+n 



TT V 

m+n 



TT 



-V m!?T,! 



^—^ ' (m -i- n — I 



^ + 1] 



fc=max(m,n) 



{m + n — k)\{k — ■m)\(k — n)\ 



(the latter formula holds for all m, n, but the formula with Lm holds for 
m < n, of course). Using the generating function [U Sect. 10.12, (17)] 

oo 

J]L(,^)(x)t- = (l-r'-^exp(-^) 

m=0 

we get also (for m < n) 



<1(0) 



TT V 




m! n-m T fn — m \ 1 d*^ 

—2— — ^ri +1^ 



m! df^ 



^0 (1 - t){n-™)/2(l + t)("-™)/2+l 



which leads (by contour integration) to a formula suitable for asymptotic 
analysis, 



<^n(0) 

(-I)*" /m! 
n 



2^ — ^ri + 1 



TT 



1 - re"^ 



dv? 



27r y (1 j-2g2i(/5^(n-m)/2+l ^^gi(py 



it holds for every r G (0, 1) and is especially useful for r = ^Jrnfn. Here are 
three asymptotic results obtained this way. 
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First, lim^^oo ^r,L+fc(0) = V(27r) for every /c = 0, 1, 2, . . . 

Second, let m,n —>■ oo, m/n c G (0, 1). If in addition m remains even, 
then w^^(O) (27r)~^c~^/^. However, if m remains odd, then u'^^(O) 
(27r)-icVl 

Third, let m be fixed and n — >• cx), then tt;^|(0) ~ Cmn^/"^ for some Cm- 
Returning to (11 ■2p we have, first, 

m,n 

= 5]cmc„<s„(0) I /(y,)e'("-'")^dy. 

for any / G i^2(0, 27r), and second. 



1.6) / fivW'^d^ 



(Sn if A; = 0, 
Q/k if (fc + 3)/12 is an integer, 
—6/k if (A; — 3) /1 2 is an integer, 
^ otherwise 

for the function / of (11.31) . The matrix elements of the operator |(i?(0) + 
E{2ti/3) + -E(47r/3)) in the basis of number states are thus calculated, 

i(m|E(0) + E(27r/3) + E(47r/3)|n) = l<s„(0) J /(^)e'('^-"^)^ d^ , 

w^^(O) being given by (II. 5p and the integral by (II. 6p . The spectrum of 
the corresponding infinite matrix determines the possible values of prob('0), 
namely, Amm < prob(V') < Amax, where Amin, Amax are the least and greatest 
elements of the spectrum. 

Restricting ourselves to a finite portion |0), |1), . . . , — 1) of the basis of 
number states, we get a finite matrix, N x N, and can compute its spectrum 
numerically. The results follow. 

TV 25 50 100 200 300 

Amax(A^) 0.6961 0.6997 0.7025 0.7045 0.7054 

Probably, the number Amax = liniAr Amax(^) is close to 0.71. But, who knows? 
It could happen that Amax = 1 and, say, 1/(1 — Xma^iN)) 3 + 0.04 In for 
large A^. 

In Section [2] it is shown that A^ax 7^ 1- 
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Here is the plot of the angular Wigner function W^^^, where ip = ^/'maxlSOO) 
is the 300-dimensional eigenvector corresponding to Ainax(300). 




Its (small) negative values are responsible for the (small) quantum violation 
of the classical bound (0.71 > 2/3). 

2 Using the Weyl transform 

We are interested in the spectral bounds of the operator 

I {E{0) + E{27r/3) + £;(47r/3)) = ^ + , 

where 

A = sgnQ(O) + sgnQ(27r/3) + sgnQ(47r/3) = 

= sgn(Q) + sgn (- + ^p) + sgn (- \q - ^p) . 

Note that {—A) is unitarily equivalent to A (since the symplectic linear trans- 
formation (g,p) ^ {—q, —p) corresponds to a unitary operator), therefore the 
spectral bounds of A are ±||>1||, and the spectral bounds of | + ^^l are 

Clearly, ||^|| < 3 (since ||sgng(t)|| = 1). 

2.1 Theorem. \\A\\ < 3. 

In order to prove the theorem, it is sufficient to prove that the spectrum 
of A does not contain 3. 

The next lemma is the first step toward this goal. 

2.2 Lemma. The number 3 is not an eigenvalue of A. 
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Proof. Assume the contrary, then 3 is also an eigenvalue of 3(| + ^A), which 
means 

{^\E{0) + E{2tt/3) + E{Att/3)\^) = 3 

for some state vector?/^, \\ip\\ = 1. Then = Ifort G {0, 27r/3, 47r/3} 

(since {ip\E(t)\tp) < 1 for these t). 

Note that [Q{0),Q(t)] = [Q, <5 cost + P sint] = isint, in particular, 
[Q(0),Q(27r/3)] = ^i. Moreover, the pair of operators (^Q(O), Q(27r/3)) 
is unitarily equivalent to the pair ((^(O), Q(7r/2)) = {Q,P) (since every 
simplectic linear transformation of the phase plane corresponds to a uni- 
tary operator). It follows that the pair of operators (E(0), E(27r/3)) = 
(^(^g(0)), 0(Q(27r/3))) is unitarily equivalent to the pair (E(0), E{tc/2)) = 

{e{Q),e{P)). 

Existence of a unit vector ip satisfying E{0)ip = = E{27T/3)tp imlies 
existence of a unit vector ipi satisfying 6{Q)ipi = ipi = 6{P)ipi. In other 
words, we get a wavefunction ipi G L2(M) concentrated on the halfline (0, cxd), 
whose Fourier transform is also concentrated on (0,oo). However, this is 
forbidden by a well-known theorem of F. and M. Riesz [5l Part One, Chapter 

1, §!]• □ 

In spite of Lemma 12.2^ the number 3 could belong to the spectrum of 
A. For example, consider the operator B = {E{0) + E{7i/2))/2 = {9{Q) + 
9{P)) /2. The number 1 is not an eigenvalue of B (recall the proof of Lemma 
12.21) . but still belongs to the spectrum of B, since {ipn\B\ipn) 1 for coherent 
states ipn such that {ipnlQli^n) +oo, {4'n\P\4'n) +00. 

In contrast, the spectrum of A is discrete, except for two accumulation 
points, —1 and 1. In order to prove this claim we prove that the operator 

— 1 is compact, moreover, belongs to the Hilbert-Schmidt class, 

trace ((A^ - 1)^) < oo . 

We do it by calculating the Weyl transform of A"^. 

In general, the Weyl transform (or Weyl symbol) of a bounded operator 
B : H —>■ H is defined as a Schwartz distribution / on such that the 
equality 

JJ fiq,p)W^{q,p)dqdp= {^p\B\^|J) 

holds for all ip of a dense set of state vectors. It is assumed that for each 
ip of this set, the Wigner function is a rapidly decreasing, infinitely 
differentiable function. It is well-known that / is uniquely determined by B. 

However, we do not need Schwartz distributions; we restrict ourselves to 
bounded measurable functions / : — > M (and the corresponding operators 
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B). Accordingly, we do not need the differentiabihty of W^. Also, we waive 
the rapid decrease of W^, demanding only J J \ W^{q,p) \ dqdp < oo. 

Clearly, the function {q,p) i— » f{aq + bp) is the Weyl transform of the 
operator f{aQ + bP), for all a, 6 G M and all bounded measurable functions 
/ : R ^ R. 

It is well-known (see [3], (3.14)) that 

(2.3) ^jj \f{q,p)\'dqdp = trace(5*5) , 

thus, an operator belongs to the Hilbert-Schmidt class if and only if its Weyl 
transform belongs to L2(M^). 

2.4 Lemma. The function {q,p) ^ (2/vr) Si{2qp) is the Weyl transform of 
the operator (sgnQ) o (sgnP). 

Here Si(x) = ^ du, and A o 5 = {AB + BA) /2. 
The (rigorous) proof grows from a non-rigorous argument shown below 
before the proof. We have 

i f°° ■ dt 
sgn Q = / e'**^ — (and the same for P) ; 

J-oo i 



2 

(sgnQ)o(sgnP) = -ij JJ ^j{<^os'j)e'^'^^^^K 

The Weyl transform of the operator e'(*'3+'*^) is the function {q,p) ^— e'*^*''"'"*^-*. 
Thus, the Weyl transform of the operator (sgnQ) o (sgnP) is the function 



7r2 



4: f f ds dt / st\ 



r2 



^ JJs>0,t>0 



I COS — ) sin tq sin sp 



s t \ 2/ 



. , ds . r dt s . 
— / — sm sp — cos -t sm qt 



rc^ Jo s Jn t 2 

ds f°^dtfl 
, — sm sp — 

Jo S Jg t ^ 

A ds . /tt / s\ Ti 
V- L T^^"^'- 4'^H^+2J + 4 



— / — sin sp / — ( - sin ( g + - ) t + - sin ( g — - ) t 



4 TT Z"^'^' ds 2 

— ■ — • sgng ■ / — sinps = — Si(2gj9) . 
TT^ 2 Jo s TT 



9 



In order to convert the calculation above into a proof we need the equality 

" ^// Tt('°4) / / dgdpe^(*^-^^^)iy^(g,p) = 

= -1 lldqdpW4q,p) II yy(cos|)e'(*^+^^). 

Unfortunately, it does not follow from Fubini's theorem, since for large s, t 
the integrand decays too slowly. (Also small s, t make a trouble, but less 
serious.) A cutoff is used below. 

Proof of Lemma 12.41 We have for all q 

/ e'*"? — = - / ^ dt = - (Si(gn) - Si(g/n)) ^ sgn q 

(as n — s> oo), and these functions are bounded in q, uniformly in n. The same 
holds for P, and we get 



(?/;|(sgnP)(sgnQ)|?/;) 



f-i)'lim hm/^f / ^e-')( [ ^e^^A 

V TtJ m^oon->oo\ \J±^,^,^^ S /VJi<|i|<n ^ / 



Therefore 



4 hm hm / - / ^fcos-)(^le^(*'3+^^)| 



{ipKsgnP) o (sgnQ)l^/') 
1 

71-2 m^oo n~.oo s Jl<\t\<n t V 

mil 71 I I 

lim lim [ — f —f cos — ) 

mil nil 



I j dgdpiy^(g,p)e^(*''+^P) 
^ lim lim // dqdpWJq,p) [ — [ ^ f cos -)e"V^P = 

71^ m-*oon^ooj J J ±^\s\<m ^ J l<\t\<n ^ ^ 2/ 

mil ri I I 

— (2i)^ lim lim / / dqdpW^{q,p) [ 

TT m— »oo n~*oo J J J ^ 



ds r dt f st^ 
, , — / — I cos — ) smgtsmps 

TT^ m^oo n~.oo 7 7 J,^^ S 7i/„ t \ 2 



However, 



dt / st\ . 
— ( cos — ) sm qt 

l/n t ^ 2 



2Jn 



is.((,.£)„)-ls.((,.£)i)H-ls.((,-£»-ls.((, 

^ -sgn -j + -sgn (^g - -j as n ^ oo 
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for all s, q, and these functions are bounded in s, g, uniformly in n. Taking 
into account that 

ff ds 

/ / dqdp\W^{q,p)\ / — | sinps| < oo 

J J Jl/m S 

we get 

(V'l(sgnP) o (sgnQ)lV') = 

= ^Ji^ jj dqdpW^{q,p) j^^ ^ sin ps- (^sgn (q+^^ + ^sgn (q-'^ 

If 2\q\ > J- then 

1^1 rrj. 



m 

^^^-smp..^-sgn(g + -j + -sgn(g-- 



TT 

2 

TT 



rmm{m,2\q\) 

sgn g) / — sin ps 

Jl/m "5 



— (sgng)(sgnp) (Si(|p| min(m, 2|g|)) — Si(|p|/m)) . 
Otherwise, if 2\q\ < the integral vanishes. We see that 

j^^ (^^^"^(^+0 +^^Sn(g-0^ ^ ^Si(2gp) as m ^ oo 

for all and these functions are bounded in uniformly in m. Taking 
into account that /J dqdp\W^{q,p)\ < oo we get 

(-01 (sgn P) o (sgn Q) IV') = jj dqdpW^{q,p) Si(2gp) . 

□ 

We return to the operators Q{t) = Qcost + Psint and recall that the 
pair of operators {Q{s), g.^^^_^^ Q(t)) is unitarily equivalent to the pair (Q, P) 
whenever sin(t — s) ^ (as was noted in the proof of Lemma [2^21 for a special 
case) . Thus, Lemma 12.41 implies the following. 

2.5 Corollary. The function 

2/2 \ 

(g, — Si —(q cos s + psin s)(q cos t + p sin t) ) 

vr V|sm()f: — s)| / 

is the Weyl transform of the operator (sgn(5(s)) o (sgnQ(t)), whenever 
sin(t - s) ^ 0. 
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2.6 Corollary. The function 



f{Q,p) = - Si 

TT 



, +-Si 



- + 



TT 



2 / 
4 — g + p\/3 —q—p\/3 
7^ 2 2 



is the Weyl transform of the operator 



-{A' - 3) = (sgng(O)) o (sgng(27r/3)) + 

+ (sgng(O)) o (sgng(47r/3)) + (sgnQ(27r/3)) o (sgn g(47r/3)) . 

2.7 Lemma. 

r I 

if{(l,p) + lydqdp < oo. 



//' 



Proof. Taking into account that / is invariant under the rotation by 7r/3 and 
the reflection {q,p) ^ {—q,p) we may integrate only over the domain < 
q\/3 < p. We divide the domain into a bounded domain (whose contribution 
is evidently finite) and two unbounded domains, one being p > 1, < q < 
p~^^^, the other p > 3, p~^/^ < q < p/VS. 
Using the inequality 



Si(x) - 



TT 



< — for < X < OO 

X 



(it holds, since 
\f{q,p) + l\ < 



smu cos a: 

X u X 



- ir W ^^)' set for < g < p/V3 



2/4 pV3- 
— Si —;=q ■ — 



1 



+ 



+ 



+ 



2/4 pV3-q pV3 + q 
'-n^'[Ts~2 2— 



) 



< 



< 



4 ^/3 



4 V3 



4 V3 



TT 4 q{pV^-q) TT 4 9(pa/3 + 9) tt 4 (p^_g)(p^ + g) 



< 



< 



v/3/ 2 2 1 

n \qp^/3 qpVS- 2/3 p^(2/3) • 3g 



< 



const 

pq 



thus 



/ dp/ dq{f{q,p) + iy < const- dp dq— 

Js Jp-y^ J3 Jp-y^ p q 



< oo . 
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Using also the inequality 

I Si(x) — Si{y) \ < \x — y\ ior x,y E 
(it holds, since |^| < 1), we get for < g < p/V^ 



\f{q,p) + i\ < 



2 / 4 pV3-q\ 2/4 pV3 + q 
— bi —;=Q ■ bi —=q ■ 



2/4 pVS-q pV^ + q 
1 bi —p^ ■ ■ 



+ 



+ 



< 



2 4 



4 V3 



4 



< q^ + 

TT V3 vr 4 (py3_g)(p73 + g) 



f 2 1 

< const ■ [q -\ — -r 



thus 



noo pp ' poo pp ' 

/ dp dg(/(g,p) + l)^ < const ■ / dp dqiq'^-l — - 

Ji Jo Ji Jo \ p 

Proof of Theorem 12.11 By (12.31) . Corollary 12.61 and Lemma E 

1 



< oo . 



□ 



trace (A^ - 1) 



2tt 



|2(/(g,p) + l)rdgdp<oo 



therefore every spectral point of A different from ±1 is an eigenvalue. By 
Lemma 12.21 the number 3 does not belong to the spectrum of A. It remains 
to use the remark after the theorem. □ 



3 Some generalizations 

3.1 Lemma. Let s, t G M satisfy s<7r<t<s + 7r. Then 

||sgnQ(0) +sgnQ(s) + sgnQ(t)|| < 3. 
Proof. Similarly to Sect. [2] we introduce the operator 

As,t = sgng(O) + sgng(s) + sgng(t) . 

Similarly to Lemma 12.21 the number 3 is not an eigenvalue of A. Thus, it 
is sufficient to prove that trace ((A^ — 1)^) < oo. Similarly to Corollary 12.61 
we calculate the Weyl transform fs^t of \{A'^^ ~ 3). It remains to prove that 
//(/s,t(g,p) + l)^dgdp< oo. 
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We observe that 



77 TT TT TT StT TT TT 

s < - <t <s + -< — <t + -<27r + s , 

2 2 2 2 2 2 2 ' 

choose tte such that 

TT TT TT TT 

ai<s--<a2<-<a3<t--<a4,<s + -< 

< as < — < < t + — < ai + 27r 

and divide the plane into six sectors by the six rays (g, p) = (r cos ak, r sin afc), 
r>0,fc = l,...,6. The three hnes q = 0, q cos s+p sin s = 0, g cos t+p sin t = 
give us six more rays, and each of these rays is contained in one (and only 
one) of the six sectors. We estimate the integral over a neighborhood of the 

ray (within the sector) similarly to / dp dq{. . . ) in the proof of Lemma 

12.71 and the integral over the rest of the sector similarly to / dp /^_i/3 dg(. . . ) 
in the proof of Lemma 12.71 □ 

3.2 Theorem. For every e > there exists 6 > such that 

||sgng(0) + sgnQ(s) + sgng(t)|| <3-5 
for all s, t satisfying s<7i — e<t — 2e<s — Se + n. 

Proof. The estimates needed in the proof of Lemma 13.11 are uniform over all 
such pairs (s,t). Thus, 

{fs,M^P) + dgdp ^0 as r oo 

uniformly in these s,t. Also, fs,t{(l,p) is continuous in s,t for any fixed q,p. 
Therefore fg^t treated as an element of L2(M^) depends continuously on s,t. 
Taking into account the isometric correspondence between Hilbert-Schmidt 
operators and their Weyl symbols we see that the Hilbert-Schmidt operator 
A^g — l depends continuously on s,t. It follows that ||y4^j — 1|| is continuous 
in s,t (the usual operator norm is meant, not the Hilbert-Schmidt norm). 
By compactness, the norm reaches its maximum at some point (sq, to) of 
the considered set of pairs (s,t). By Lemma IXTl HAsQ^fgll < 3. Therefore 
Pso to-M\ < 8- We take S > such that (3 - 5)^ - 1 > \\Al^ to - M\ and 
(3 -5)2 - 1 > 1. Then WA^^ - 1|| < (3 - 5)^ - 1 implies \\As,t\\ < 3 - 5. □ 
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